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^ . 

^' A local turbulence model is developed to study energy cascades in the inter- 

tyj ' stellar medium (ISM) based on self-consistent two-dimensional fluid simulations. 

The model describes a partially ionized magnetofiuid interstellar medium (ISM) 
J>-,' that couples a neutral hydrogen fluid with a plasma primarily through charge ex- 

change interactions. Charge exchange interactions are ubiquitous in warm ISM 
plasma and the strength of the interaction depends largely on the relative speed 
between the plasma and the neutral fluid. Unlike small length-scale linear coUi- 
sional dissipation in a single fluid, charge exchange processes introduce channels 
cn . that can be effective on a variety of length-scales that depend on the neutral and 

^r^ ■ plasma densities, temperature, relative velocities, charge exchange cross section 

[--. . and the characteristic length scales. We find, from scaling arguments and non- 

linear coupled fluid simulations, that charge exchange interactions modify spec- 
O ■ tral transfer associated with large scale energy containing eddies. Consequently, 

the warm ISM turbulent cascade rate prolongs spectral transfer amongst iner- 
S^ . tial range turbulent modes. Turbulent spectra associated with the neutral and 

plasma ISM fluids are therefore steeper than those predicted by Kolmogorov's 
phenomenology. 
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Introduction 



The local interstellar medium (LISM) surrounding the heliosphere is warm (~ TOOC-ft') 
and consequently par tially ionized. Althoug h the ioni zation fraction of the LISM is not co n- 
clusively estabhshed ( ISlavin fc Frisch Il2006l )[see also (JMuUer et al II2006I : IZank et al II2006I )]. 
it is possible that the plasma number density is as small as 0.05cm~^ and the neutral atomic 
H number density as much as 0.16cm~^. Indeed, much of the interstellar medium is warm 
and partially ionized. In the LISM, the low density plasma and neutral H gas are coupled 
primarily through the process of charge exchange. On sufficiently large temporal and spatial 
scales, a partially ionized plasma is typically regarded as equilibrated; this is the case for 
the LISM. H owever, the region between the bow shock of our heliosphere and the heliopause 
( Zank II1999I ) is not equilibrated because the charge exchange mean free path (mfp) and 
the size of the region are comparable. The neutral H distribution in the region is com plex 



( iPauls et. al I Il995l : I Zank fc Pauls I Il996l : IZank et al. I Il996l : iHeerikhuisen et al I I2OO6I ) pos- 



sessing an approxima tely Maxwel 
the in ner heliosheath (IZank et al. 



ian c ore that is broadened by hot H atoms produced in 



19961 ) and fast neutral H produced in the supersonic solar 



wind ( Zank et al. Ill996l : IZank Ill999l ). Similarly, the atmospheres of many stars, especially 



cooler stars, that are embedded in a partially ionized interstellar medium will have outer 
astrosheaths that are partially ionized. Several examples of such parti ally ionized astrop 



sheres have been observed using Lyman-a absorption measurements (JGayley et al 
Woodet. allboooh . 



1996 



Turbulent fluctuations in a partially ionized LISM are not only potentially important 
in the context of the global heliospheric ISM interaction, but are also instrumental to our 
understanding of many astrophysical phenomena including the energization and transport of 
cosmic rays, gamma-ray bursts, ISM density spectra etc. Besides the ISM, a partially ionized 
plasma environment represent an important component of the Satur nion magn etosphere, es- 
pecially that part infl uenced by mass loss f rom the moon Bucephalus (jlp 1119971 ). in Neptune's 
magnetosphere e.g., (JHill fc Dessler Ill990l ). and also in the astrosphere of stars embedded in 
a cloud of neutral H or even some stellar winds. For magnetosphere, the partially ionized 
plasma sometimes comprises a mixture of water group and O neutrals and pickup ions and 
solar wind pla sma. The partially ion ized plasma environment is also relevant to the study of 
neutron stars (jPotekhin et al II2005I ) where the plasma and strong magnetic field interaction 
still remains an unresolved issue. 

The interaction of a neutral gas and plasma can be dated back to the seminal work 
by Kulsrud & Pearce (1969) in the context of cosmic ray propagation, where it was shown 
that neutral component damps Alfven waves. Neutrals interacting with plasma via a rela- 
tive drag process results in ambipolar diffusion. Ambipolar diffusion plays a crucial role in 
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the dynamical evolution of the near solar atmosphere, interstellar medium, and molecular 
clouds and star formation. For instance, Oishi & Mac Low (2006) investigated the inability 
of ambipolar diffusion to set a characteristic mass scale in molecular clouds and find that 
substantial structure persists below the ambipolar diffusion scale because of the propagation 
of compressive slow mode MHD waves at smaller scales. It was further shown that the spec- 
tral behaviour is not influenced by the ambipolar diffusive process unlike viscous damping. 
Leake et al (2005) showed that the lower chromosphere contains neutral atoms, the existence 
of which greatly increases the efficiency of wave damping due to coUisional friction momen- 
tum transfer. They find that Alfven waves with frequencies above 0.6Hz are completely 
damped and frequencies below 0.01 Hz are unaffected. Khodachenko et al (2006) under- 
took a quantitative comparative study of the efficiency of the role of (ion-neutral) collisional 
friction, viscous and thermal conductivity mechanisms in damping MHD waves in different 
parts of the solar atmosphere. It was pointed out by the authors that a correct description of 
MHD wave damping requires the consideration of all energy dissipation mechanisms through 
the inclusion of the appropriate terms in the generalized Ohms law, the momentum, energy 
and induction equations. Molecular clouds are known to be supported, at least in part, by 
magnetic fields. The removal of magnetic fields thus represents an important component of 
the star formation process. In the most studied scenario, field removal occurs through the 
action of ambipolar diffusion, wherein magnetic fields are tied to the ionized component, 
which drifts relative to the more dominant neutral component of the gas (Mestel & Spitzer 
1956, Mouschovias 1976, Nakano 1979, Shu 1983, Nakano 1984, Lizano & Shu 1989, Basu & 
Mouschovias 1994). The role of magnetic fields and ion-neutral friction in regulating gravita- 
tionally driven fragmentation of molecular clouds was studied by Kudoh et al (2007). Mestel 
& Spitzer (1956) pointed out that even if clouds are magnetically supported, ambipolar dif- 
fusion (resulting from ion-neutral drag) will cause the support to be lost and stars to form. 
Li & Nakamura (2004) and Nakamura & Li (2005) have studied turbulent fragmentation 
processes for a magnetized sheet including the effect of ion-neutral friction. They find that 
a mildly subcritical cloud can undergo locally rapid ambipolar diffusion and form multiple 
fragments because of an initial large scale highly supersonic compression wave. Padoan et al 
(2000) calculated frictional heating by ion-neutral (or ambipolar) drift in turbulent magne- 
tized molecular clouds and showed that the ambipolar heating rate per unit volume depends 
on field strength for constant rms Mach number of the flow, and on the Alfvenic Mach num- 
ber. Furthermore, ion-neutral friction, has long been thought to be an important energy 
dissipation mechanism in molecular clouds, and therefore a significant heating mechanism 
for molecular cloud gas (Scalo 1977, Goldsmith & Langer 1978, Zweibel & Josafatsson 1983, 
Elmegreen 1985). These are only a few, amongst numerous other, studies that illustrate 
the importance of a neutral gas component in the dynamics of partially ionized magnetized 
astrophysical plasmas. 
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On even smaller scales, the edge region of a tokamak (a donut shaped toroidal experi- 
mental device designed to achieve thermonuclear fusion reaction in an extremely hot plasma) 
is partially ionized. In a tokmak, the neutral particles result from effects such as gas puff- 
ing, impurity injection, recombination, charge exchange, and possibly neutral beam injection 
processes. The presence of neutrals can potentially alter the dynamics of zonal flows and 
cross-field diffusivity. For instance, several experiments have demonstrated that the tran- 
sition from L-mode (low confinement) to H-mode (high con finement) can be significantly 



affected by neutral atoms in the edge of a tokamak plasma (jSingh et al II2004I ). The edge 
region predominantly consists of neutral species resulting from recycling from the wall, and 
from the limiter and/or divertor plate. These neutrals are present in significant numbers (e.g. 
10^° — 10^^ cm~^, which is about 5-8% of the plasma density near the edge) and affects the 
poloidal momentum balance and the L-H transition threshold. Thus they a play crucial role 
in regulating transport processes in tokamak fusion plasmas. The neutral gas can interact 
with the magnetized plasma via charge exchange and other effects. Furthermore, Doublet 
III-D tokamak (DIII-D) data shows that there is a significant correlation between the power 
threshol d for the L-H transitio n and the poloidally averaged neutral density at the 95% flux 



surface ( JGroebner et al Il200ll ). To understand the role of neutral gas on H-mode physics. 



asymme tric gas puffing experimen ts have also been performed in COMPASS-D and MAST 



devices f IValovie et al Il200ll . 120021 ). 



In all these environments, turbulence is ubiquitous. Turbulence involves the nonlinear 
coupling and transfer of energy across a multitude of spatial and temporal scales. The 
coupling of plasma to a neutral gas via charge exchange, especially in a non-equilibrated 
environment, introduces both a length scale distinct from the usual collisional mfp and 
alternate channels and mode couplings for the transfer of energy. Certainly, in the local 
ISM, the charge exchange cross-section of e.g. neutral hydrogen atoms is larger than the 
collisional cross-section (e.g. Zank 1999), implying that the former governs much of the 
small-scale physics. Fluctuations with length scales that exceed the charge exchange mean 
free path ice will obviously be effected by the neutral gas, such as the damping of linear 
modes (Shaikh & Zank 2006), but it is less clear how turbulent plasma motion and nonlinear 
couplings will be mediated. For plasma fluctuations on scales smaller than i^e, the role of 
neutral gas is even less clear. 

Further, it is worth mentioning that the charge exchange process in its simpler (and 
leading order) form can be treated like a friction or viscous drag term in the fluid momen- 
tum equation, describing the relative difference in the ion and neutral fluid velocities. The 
drag imparted in this manner by a collision between ion and neutral also causes ambipolar 
diffusion, a mechanism used to describe the Alfven wave damping by cosmic rays (Kulsrud 
& Pearce 1969) and also discussed by Oishi & Mac Low (2006) in the context of molecular 
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clouds. The thermodynamical properties of the latter are significantly different from that 
found in the heliosheath region. Interestingly, unlike linear damping, ambipolar diffusion 
does not terminate an isothermal MHD turbulence cascade (Oishi & Mac Low 2006). In- 
stead, it damps some linear MHD waves. In the context of MHD turbulence, this issue had 
also been addressed by Cho & Lazarian (2003) who discussed the viscosity-damped effects 
on inertial subrange associated with magnetic field fluctuations in compressible (isothermal) 
MHD turbulence. Ion-neutral friction nonetheless describes a leading order interaction be- 
tween the ion and neutral fluid species where thermal speeds associated with both the fluids 
are of little significance or ignorable and plasma phenomena are describable predominantly 
by isothermal processes. This description is not appropriate to describe the outer heliosheath 
plasma where proton as well as neutral hydrogen temperatures can be as high as 20, 000 °K. 
The thermal speeds associated with the heliosheath protons and neutrals corresponding to 
such a high temperature cannot therefore be ignored for the following reasons. Firstly, since 
the thermal speed describes a significant fraction of the energy and is proportional to the 
relative speeds in the momentum and energy charge exchange sources, it should not be 
neglected for plasmas that are not cold. It can potentially influence the momentum and 
energy sources during the course of evolution. This means that the interaction between the 
heliosheath protons and neutral atoms is not mediated only by the corresponding bulk fluid 
speeds, but it also depends critically on their thermal speeds. Secondly, for elastic and charge 
exchange collisions, the initial particle distribution can be well described by Maxwellians even 
with different streaming velocities and temperatures. Higher order effects are introduced in 
the dynamics because the distributions are not really Maxwellians due to prior collisions as 
pointed out by McNutt et al. (1998). The latter introduces Navier-Stokes-like corrections 
(viscosity and thermal conductivity like terms) to the neutral component of the equation. 
A consistent treatment of such processes, describing the higher order corrections, in the 
context of the heliosheath requires that the Maxwellians be multiplied by transfer integrals, 
as done in Pauls et al (1995). Thirdly, plasma and neutral fluids in the heliosheath region 
are compressible enough to possess a finite component of thermal energy associated with 
a non-adiabatic exchange of energy amongst the local fluctuations. Thus in the context 
of supersonic heliosheath turbulent fluctuations, the characteristic speed of plasma-neutral 
coupled fluid needs to be treated on equal footing with the thermal speed. Consequently, 
thermal corrections appear as higher order terms along with the leading order ambipolar-like 
diffusive forces in the charge exchange expressions as described in section 3. This will be 
evident from the full energy equation and the subsequent charge exchange energy transfer 
functions employed in our model. 

Multiple processes can couple the interstellar medium plasma and neutral gas. Besides 
charge exchange scattering due to an induced dipole moment in the neutral atom, depending 
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on the temperature (Banks, 1966), is possible. Interestingly enough, variations in cross- 
sections can also be expected for different ions in the same neutral gas at higher temperatures 
(Banks, 1966). Further processes, especially in the context of heliospheric processes, are 
photoionization, electron impact ionization, recombination, and H-H"*" Coulomb collisions 
(see e.g. Zank, 1999). However, cross-section associated with most of them is small compared 
to that of charge exchange between hydrogen atom and ions in the heliosheath region. We 
therefore do not incorporate them in our current model. 

In view of these comments, we have developed a self-consistent plasma-neut ral ISM tur- 



bulence simulation model based on fl uid modeling of neutral and plasma fluids (jZank Ill999 



Pauls et. al Ill995l : IZank et al. Ill996l ). Our model treats the plasma and neutral species as 
distinguishable fluids and couples them through charge exchange. The model essentially 
employs a set of magnetohydrodynamic equations that describe the turbulent plasma mo- 
tions, whereas the neutrals are described by hydrodynamic equations. A fluid approach for 
the local interstellar medium is justifiable since the charge exchange and collisional mean 
free path ensures approximate thermodynamical equipartition between plasma and neutral 
species with the result that both plasma and neutrals can be described by a Maxwellian dis- 
tribution. Besides the local ISM, the outer heliosheath can also be descri b ed re asonably well 
by a fluid model of this kind, as was done in the model of iPauls et. al I (119951 ). For a non- 
eq uilibrated descripti o n, a more elabo r ate ra ulti-fluid model can be developed, as was done 
by lZank et al. I (119961 ): IZank fc Pauls I (119961 ). and this could be applied to straightforwardly 
to both the outer heliosheath and very local ISM. This would be more important to the outer 
heliosheath since a multi-fluid model would then include hot inner heliosheath neutral atoms 



1996 



deposited through secondary cha rge exchange in the outer heliosheath (IZank et al. 

Zank fc Pauls Ill996l : IZank Ill999l ). The schematic in Fig. ([1]) depicts the simulation regions 



in which we are particularly interested i.e. the outer heliosheath and the local interstellar 
medium. The importance of ISM turbulence is clearly evident from this figure in the context 
of cosmic ray modulation, large-scale heliospheric structure and particle acceleration. 

In this paper, we focus on the multi-scale evolution of a partially ionized plasma, com- 
prising multiple fluids, where each fluid evolves under the influence of the other through a 
complex interaction process. We investigate the energy cascade between LISM turbulent- 
fluctuations in a partially ionized plasma where the plasma and neutral gas are coupled 
through charge exchange. Here, we, present a numerical simulation of ISM turbulence that 
self-consistently evolves both the plasma and neutral fluids when coupled by charge ex- 
change. We restrict our attention to a neutral hydrogen (H) gas, so that plasma and neutral 
densities are conserved by charge exchange interactions. Hence the corresponding continuity 
equations do not include charge exchange terms. One of the points that emerges from our 
coupled multi fluid LISM simulations is that neutrals enhance turbulent cascade rates and 
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l ead to much steep er spectra compared to that predic ted by the usual Kolm ogorov theory 
( Kolmo gorovl 1 1 94 ll ) for pure MHD or hydrodynamics ( Shaikh &: Zank 1120071 ). In Section 2, 



we discuss the equations of a coupled plasma-neutral ISM turbulence model, their validity, 
the underlying assumptions and the normalizations. Section 3 deals with a quantitative 
derivation of charge exchange source terms, associated with the plasma and neutral fluids, 
based on the time-dependent Boltzmann equation. Section 4 describes the results of our 
nonlinear, coupled, self-consistent ISM fluid simulations. It is shown, in section 5, that the 
coupling of plasma and neutral fluids leads to an enhanced spectral transfer in the inertial 
range, which results in a steeper ISM turbulent spectrum than predicted by the Kolmogorov 
theory. Finally, conclusions are presented in Section 6. 



2. Model Equations for a partially ionized LISM 

The assumptions that are intrinsic to our model of turbulence in a partially ionized 
ISM are the following, (i) Fluctuations in the plasma and neutral fluids are isotropic, ho- 
mogeneous, thermally equilibrated and turbulent, and (ii) Neither a mean magnetic field 
nor velocity flows are present initially. Local mean flows may subsequently be generated 
by self-consistently excited nonlinear instabilities, (iii) The characteristic turbulent corre- 
lation length-scales (Ac ~ 1//Cc) are typically bigger than charge-exchange mean free path 
lengths (Ace ~ 1/^ce) in the ISM flow s , i.e A^ ^ Ace or kce/kc 3> 1. The latter inequality 



is also consistent with iFlorinski et al I (120031 . |2005| ). Nevertheless, they are large enough to 



treat any localized shocks as smooth discontinuities. In other words, the characteristic shock 
length-scales are small compared to the ISM turbulent fluctuation length-scales, and finally 
(iv) boundary conditions are periodic, essentially a box of interstellar plasma, which is a 
natural and most appropriate choice for modeling turbulence in the local ISM. 

While most of the above assumptions are appropriate to realistic ISM turbulent flows, 
so allowing us to use MHD and hydrodynamic descriptions for the plasma and the neutral 
components respectively, the latter has been criticized in the context of heliospheric flows 
within the heliopause owing to the large charge-exchange mean free path (about the size or 
even bigger than the extent of the heliosphere). The use of a fluid description for neutrals in 
the inner heliosheath is therefore inappropriate. In the LISM and outer heliosheath, we are 
nonetheless not restricted by a large mfp because the plasma and neutral fluid remain close 
to thermal equilibirium and behave as Maxwellian fluids. Our model simulates the local 
ISM and outer heliosheath. The fluid model describing nonlinear turbulent processes in the 
interstellar medium, in the presence of charge exchange, can be cast into plasma density {pp)., 
velocity (Up), magnetic field (B), pressure {Pp) components according to the conservative 



form 



where, 



at 



+ V • Q„ = Q, 



■p =i^p,ni 



Pp 

PpUp 
B 



,Qp 



^v.n 



and 



PpUp 
PpU.U, + :^ + f - BB 

UpB - BUp ' ^^ 

epUp - B(Up ■ B) 



1 rr2 Pp B^ 





Qa/ (Up, V „, Pp, p„, r„, Tp) 



Q_E(Up, V „, Pp, p„, T„, Tp) 



The above set of plasma equations is supplemented by V • B = and is coupled self- 
consistently to the ISM neutral density (p„), velocity (V^) and pressure (P„) through a set 
of hydrodynamic fluid equations. 



at 



+ V ■ Q„ = Q. 



n,p} 



(2) 



where. 



Pn 
Pn * r. 



^r 



Pn * 71 
n V V + -^ 



=£-71,^ 





'^MI^ "^ n; '-'pj Pp; Pn; -^ 71; -'pj 
^1 E\^ ni L'p; Pp; P7I; -'71; -tpj 



1 P 

en = -PnK + -. 

/ 7 — 1 

Equations ([I]) to ([2]) form an entirely self-consistent description of the coupled ISM 
plasma-neutral turbulent fluid. Several points are worth noting. The charge-exchange mo- 
mentum sources in the plasma and the neutral fluids, i.e. Eqs. ([1]) and ([2]), are described 
respectively by terms QA/(Up, V„, pp, pn, Tn, Tp) and Qa/ (V„, Up, pp, pn, Tn, Tp). A swapping 
of the plasma and the neutral fluid velocities in this representation corresponds, for instance, 
to momentum changes (i.e. gain or loss) in the plasma fluid as a result of charge exchange 
with the ISM neutral atoms (i.e. Qm(Up, Vn,Pp,Pn,Tn,Tp) in Eq. ([1])). Similarly, mo- 
mentum change in the neutral fluid by virtue of charge exchange with the plasma ions is 
indicated by Qm(V„, Up,Pp,p„,T„,Tp) in Eq. ([2]). In the absence of charge exchange inter- 
actions, the plasma and the neutral fluid are de-coupled trivially and behave as ideal fluids. 
While the charge-exchange interactions modify the momentum and the energy of plasma and 
the neutral fluids, they conserve density in both the fluids (since we neglect photoionization 
and recombination). Nonetheless, the volume integrated energy and the density of the entire 
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coupled system will remain conserved in a statistical manner. The conservation processes 
can however be altered dramatically in the presence of any external forces. These can in- 
clude large-scale random driving of turbulence due to any external forces or instabilities, 
supernova explosions, stellar winds, etc. Finally, the magnetic field evolution is governed 
by the usual induction equation, i.e. Eq. ([1]), that obeys the frozen-in- field theorem unless 
some nonlinear dissipative mechanism introduces small-scale damping. 

The underlying ISM turbulence model can be non-dimensionalized straightforwardly 
using a typical ISM scale-length {£q), density (po) and velocity (vq). The normalized plasma 
density, velocity, energy and the magnetic field are respectively; pp = Pp/po, Up = Up/vo, Pp = 
Pp/poVQ, B = B/fo^po- The corresponding neutral fluid quantities are p„ = Pn/ Po, U„ = 
\Jn/vQ, Pn = Pu/poVq- The momentum and the energy charge-exchange terms, in the nor- 
malized form, are respectively Qm = Qm£o/poVQ,Qe = Qc^o/poVq- The non-dimensional 
temporal and spatial length-scales are i = tvo/io, x = x/£o- Note that we have removed bars 
from the set of normalized coupled ISM model equations ([I]) & ([2]). The charge-exchange 
cross-section parameter (a), which does not appear directly in the above set of equations (see 
the subsequent section for more detail), is normalized as a = uoioa, where the factor rioio 
has dimension of (area)~^. By defining no,i'o through ace = ^/noio = k~^, we see that there 
exists a critical charge exchange wavenumber (fcce) associated with the coupled ISM plasma- 
neutral turbulent system. For a characteristic density, this corresponds physically to an area 
defined by the charge exchange mode being equal to (mfp)^. Thus the larger the area, the 
higher is the probability of charge exchange between plasma ions and neutral atoms, as illus- 
trated in Fig. (2). To be precise, 'fcce' is a (lengthscale)"^ that typically helps us determine 
whether or not a particular turbulent fluctuation length scale undergoes charge exchange in 
our model. It is strictly in that sense, we refer to it as a critical wavenumber. Clearly, there 
exist three different regimes depending on whether (i) k < kce, (ii) k ~ kce or (iii) k > kce- 
In the heliosheath, the probability that charge exchange can directly modify those modes 
satisfying k < kce is high compared to modes satisfying k > kce- Since the charge exchange 
length-scales are much smaller than the ISM turbulent correlation scales, this further allows 
many charge exchange interactions amongst the nonlinear turbulent ISM modes before they 
cascade energy in one unit convective time. This is illustrated schematically in Fig. (2). 

It is interesting to compare our model with other work. Of particular relevance is the 
classic work by Kulsrud and Pearce (1969) who investigated the interaction of galactic comic 
rays and Alfien waves in the interstellar medium and also considered ion-neutral collision as 
a damping mechanism for Alfven waves. Ion-neutral collisions are also reported by Mclvor 
(1977) to be a dominant mechanism for damping interstellar medium turbulence. The damp- 
ing (due to ion-neutral collisions) rates for Alfven waves estimated by Kulsrud and Pearce 
(1969) are based on uq ~ noVthO', where z/o,no,fth and a are respectively the ion-neutral 
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collision frequency, neutral density, thermal speed and collision cross-section. The damping 
wavenumber associated with this frequency can be estimated to be /c ~ nQVth<^ /Va- The 
latter determines essentially the dissipation wavenumber associated with the damping of 
small scale Alfven waves in ISM turbulence (Kulsrud and Pearce 1969). The turbulent cas- 
cade is further expected to be terminated beyond this wavenumber. It is clear from this 
expression that the higher the thermal speed, the larger the dissipation wavenumber. Corre- 
spondingly in real space, dissipation will be concentrated on smaller length-scales for higher 
thermal speeds. Furthermore, the thermal speed is also proportional to the temperature. 
Thus, given the heliosheath parameters, vthu associated with heliosheath neutral atoms will 
be at least two orders of magnitude larger compared with that results from Eq. (47) in 
Kulsrud and Pearce (1969). This will correspond to a larger dissipative wavenumber in the 
heliosheath region and essentially means that charge exchange interactions are restricted 
to only those eddies whose length-scales are comparable to the dissipative length-scales of 
turbulence. It is not clear if this scenario holds in inner /outer heliosheath turbulence where 
charge exchange interactions, in principle, can modify any length scale in the inertial range. 
Part of the reason resides with the fact that charge exchange modifications are proportional 
to the relative speed between ion and neutral fluids and corresponding densities which can 
obviously alter the nonlinear cascades by influencing the combined momenta and energies 
of both the fluids in a self-consistent and subtle manner. Additionally, Kulsrud and Pearce 
(1969), as does Mclvor (1977), assumed a flxed neutral density while computing the damp- 
ing rates associated with Alfven waves in ISM turbulence. While a flxed neutral density 
can modify the linear growth rates of the underlying waves, it does not itself evolve and 
hence no back reaction on Alfven waves and vice versa can be expected. By contrast, our 
model self-consistently evolves neutral and plasma and describes the mutual feedback of one 
species on the other. This not only alters the linear interaction process, but it also influences 
the nonlinear coupling in a subtle manner not readily describable by means of any linear 
analytic theory. A self-consistent coupling is absolutely essential to understand nonlinear 
turbulent cascades in the inner /outer heliosheath. It is to be noted further that we do not 
attempt to investigate damping of heliosheath turbulence by Alfven waves that experience 
ion-neutral collisions, unlike Mclvor (1977). While Alfven waves are intrinsically present in 
our model, turbulent cascades are least affected by them along the mean magnetic fleld (She- 
balin & Montgomery 1983). Alfven waves can be crucial in providing a local anisotropy in 
the spectral transfer due to a mean or local fleld though. However, our prime focus here is to 
develop an understanding of heliosheath turbulence where plasma and neutral fluids evolve 
through a mutual interaction mediated by charge exchange. As evident from the complexity 
associated with these interactions, nonlinear turbulent cascades are likely to be affected at 
any length scale in the heliosheath inertial range, unlike a linear dissipative processes. 
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In the subsequent section, we derive an exact quantitative form of the sources due to 
charge exchange in ISM turbulence. 



3. Charge Exchange Sources 

The charge exchange terms can be obtained from the Bohzmann transport equation 
that describes the evolution of a neutral distribution function /„ = /(x, v^-, Vy, v^, t) in a six- 
dimensional phase space defin ed respectively by p osition and velocity vectors (x, t>^,t>j^,t>2) 



at each time t. Here we follow iPauls et. al I (119951 ) in computing the charge exchange terms 



from various moments of the Boltzmann equation. The Boltzmann equation for the neutral 
distribution contains a source term proportional to the proton distribution function fp and 
a loss term proportional to the neutral distribution function /„, 

-^ + V„ • Vfn H Vv„/n = /p(x, V^, Vy, V^, t) / /„(x, V^, Vy, V^, t)|f„ - Up\ae^{Vrel)d^Vn 

-/„(x. „.. ... v., t) j /,(x. „.. ... „.. t)K - .>.,(.„,)<*'«. (3) 

Here /p, Up represent respectively the proton distribution function and velocity. Uex is the 
charge exchange cross-section (between neutrals and plasma protons), m is the mass of par- 
ticle, and F represents forces acting on the fluid. The charge exchange parameter has a 
logarithmically weak dependence on the relative speed [ vrei = l^p "7 '^rX) of the neutrals and 
the protons through o"ex = [(2.1— 0.092 ln(t>re;))10~^cm]^ ( iFite et al Ill962l ). This cross-section 



is valid as long as energy does not exceed \eV , which usually is the case in the inner /outer 
heliosphere. Beyond \eV energy, this cross-section yields a higher neutral density. This issue 
is not applicable to our model and hence we will not consider it here. The density, momen- 
tum, and energy of the thermally equilibrated Maxwellian proton and neutral fluids can be 
computed from Eq. ([3]) by using the zeroth, first and second moments / f^d^^, f ^f^d^C, 
and J m^'^/2f^d^^ respectively, where ^ = Up or f„. Since charge exchange conserves the 
density of the proton and neutral fluids, there are no sources in the corresponding continuity 
equations. We, therefore, need not compute the zeroth moment of the distribution function. 
Computing directly the first moment from Eq. ([3]), we obtain the neutral fluid momentum 
equation as given by Eq. ([2]) . The entire rhs of Eq. ([3]) can now be replaced by a momentum 
transfer function QM{vn,Up) which reads 

QhliVn, Up) = p,{Up, Vn) - /i(f„, Up), (4) 

where Qm and /2, the transfer integral, are vector quantities. The transfer integrals describe 
the charge exchange transfer of momentum from proton to neutral fluid and vice versa. The 
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first term on the rhs of Eq. (4) can be expressed by 

where, 

/3(.,.„„) = //„(x.,...,.„..*)|...-.,|.„(.„,)A„. 

Considering a MaxweUian distribution for the neutral atoms, we simphfy (3{up, Vn) as follows. 



p{up, Vn) = a^^nnVr^ 



4 {Vn- UpY 



+ 



vr 






Note that the above expression emerges directly from a straightforward integration of sources 
in the rhs of the Boltzmann Eq. (3). To obtain the expression for the momentum transferred 
from proton to neutral (or vice versa), we need to take a second moment of the fl{up,Vn) 
expression. This is shown in the following, 

fi{Up, Vn) = m\nh{Up, Vn) + m{Up - V„)/i(Mp, Vn). 

where Jq and Ji are transfer integrals that can be written as follows, 

hiUr^^Vr,) = I fp{:>C,Vx,Vy,Vz,t)(3{Up,Vn) d^Up] 



lO{Up,Vn) 



h{Up,Vn) 



Vnfp(;>^,Vcc,Vy,Vz,t)P{Up,Vn) d^Up. 



Assuming a MaxweUian distribution for plasma protons and using (3{up, f„) from the above 
expression, we can straightforwardly evaluate the transfer integrals Jq and Ji (see Pauls et 
al (1995) for details). We further write the complete form of the first term on the rhs of Eq. 
(4) as follows. 



IJ,{Up,Vn) = ma^yjipHn 



u* 






ln,-v„ 



"^p t'^n 



In a similar manner, we can evaluate the second term on the rhs of Eq. (4), which yields the 
following form. 



where 5Vu v 



u* = u: 



fJ.{Vn,Up) = maeyjipUn 
4 Uvl + Af/2) + ^-^Vl 



ULu^Up - (Vn - Up 



V^ 



OVv„,Up 



4 '^T^ 



1/2 



97rT/2 ll/2 



U* 



±VS + ^VS + Af/2 



,5K„,..= [4(il^^„+At/^)+fV-J^^^and 

ll/2 I II 1 

, Am = Up — Vn (iPauls et. al 1 119951 ) . On 



substituting these expressions in the momentum transfer function, we obtain 



QM{Vn,Up) = maeyJlpUniyn " % 



U* + 



V, 



Tn 



V^ 



SVu„,v„ sv^. 



(5) 
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Equation ([5]) together with Eq. ([2]) yields the momentum equation for the neutral gas. 
Swapping the plasma and neutral fluid velocities yields the corresponding source term for the 
proton fluid momentum equation. The gain or the loss in neutral or proton fluid momentum 
depends upon the charge exchange sources, which depend upon the relative speeds between 
neutrals and the protons. The thermal speeds of proton and neutral gas in Eq. ([5]) are given 
respectively by V^ = 2KBTp/m (the factor 2 arises because of thermal equilibration in that 
it is assumed that the temperature of the plasma electrons and protons are nearly identical 
so that Tp = Te + Tproton — 2Tp) and V^^ = KBTn/m. The corresponding temperatures 
are related to the pressures by Pp = 2npKBTp and P„ = UnKsTn, where nn,Tn,np,Tp are 
respectively the neutral and plasma density and the temperature, and Kb is the Boltzmann 
constant. 

The moment, J 'mC,'^/2f^d^^, of the Boltzmann Eq. ([3]) yields an energy equation for 
the neutral fluid whose rhs contains the charge exchange energy transfer function 

QE{Vn, Up) = r]{Up, Vn) - vi^n, Up), 

where r]{up,Vn), rj{vn,Up) are the energy transfer (from neutral to proton and vice versa) 
rates. These functions can be computed as follows: 

13 V"^ 

r]{up, Vn) = -mVnCr^^npnJJ* + -mV:^^aexnpri„AKp,,;„ - ma^^UprinVn ■ (Up - v„) '^" 

^4 Vup,Vn 

and 

1 3 ^2 

r]{vn, Up) = -mV^anpHnU* + -mV^^(Tcxr^p^nAK„,nj, - ma^^npUnUp ■ (v„ - Up) 



The total energy transfer from neutral to proton fluid, due to charge exchange, can then be 
written as, 

-'- -r*/'T/2 tt2\ "J 



QE{vmUp) = -ma^^npnrJJ*{V^ - Up) + -mcrexnpn„(l^y^AKp,i,„ - V^^/Wy^^up 



-mae^npHn 



Up - v„)— — ^^^ Up ■ (v„ - Up) 



' n V P 



6V "^ ' '^'SV 



(6) 



with AVu V = \-V^ + ^V^ + At/2]i/2 and AV, , = [^V^ + ^V^ + At/^lVs. a similar 

"pi"n Itt J-p 9tt J-n J Vn,u,p I TT -^ n 9tt J- p J 

expression for the energy transfer charge exchange source term of plasma energy in Eq. ([T]) 
can be obtained by exchanging the plasma and neutral fluid velocities. In the normalized 
momentum and energy charge exchange source terms, the factor maex in Eqs. ([5]) & ([6]) is 
simply replaced by a. 
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4. Nonlinear Simulation of Energy Cascades 

It is clear from previous sections that the interaction of the neutrals with a plasma 
introduces characteristic length and time scales that can separate characteristic time and 
length scales in ISM turbulence. For instance, the plasma-neutral fluid coupling in the ISM 
fluctuations introduces charge exchange modes, kce, that are distinctively different from the 
characteristic turbulent mode k. Typically, characteristic turbulent length-scales (Ac ~ l/A;) 
are much bigger than the charge e xchange length-scales A^p, ~ I /kce, since for example, Ac ~ 
(100 - lOOOU.. in the ISM cloud Jpiorinski et al Iboosl . boosh. This essentially corresponds 



to kce/k > 1 (IZanklll999l : iFlorinski et al 



2003 



20051 ) ■ It should be noted that the charge 



exchange sources in the neutral and the plasma fluids do not merely damp the smallest scale 
fluctuations, such as is the case with linear coUisional dis sipation in turbulent fluid flows, but 



they can alter the nonlinear cascade rate dramatically (jShaikh et al I l2006l : IShaikh &: Zank 
20071 ). Because of its complex nonlinear character (see Eqs. dS]) & ([6])), charge exchange can 



be effective on a variety of turbulent length-scales unlike linear dissipation which is operative 
on small length-scales only. The effect and influence of the charge exchange source terms on 
the turbulent inertial range of coupled neutral and plasma fluids is therefore a subtle issue 
and is not restricted only to the damping of small scales in partially ionized, magnetized 
ISM turbulence. With this in mind, we carry out two different sets of simulations of Eqs. 
([T]) & ([2]), to discern the effects of charge exchange on the ISM plasma-neutral system. In 
the flrst set of simulations, Eqs. ([I]) & ([2]) are decoupled (i.e. no charge exchange coupling is 
assumed) and each evolves as an ideal fluid, whereas in the second set the coupling is included 
self-consistently The former provides a reference simulation against which to investigate the 
coupled system. 

A two-dimensional (2D) nonlinear fluid code was developed to numerically integrate Eqs. 
([1]) to ([2]). The 2D simulations are not only computationally less expensive (compared to a 
fully 3D calculation), but they offer signiflcantly higher resolution (to compute inertial range 
turbulence spectra) even on moderately-sized clusters like our local Beowulf. The spatial 
discretization in our code uses a discrete Fourier representation of turbulent fluctuations 
based on a pseudospectral method, while we use a Runge Kutta 4 method for the temporal 
integration. All the fluctuations are initialized isotropically (no mean flelds are assumed) 
with random phases and amplitudes in Fourier space. This algorithm ensures conservation 
of total energy and mean fluid density per unit time in the absence of charge exchange 
and external random forcing. Additionally, V ■ B = is satisfled at each time step. Our 
code is massively parallelized using Message Passing Interface (MPI) libraries to facilitate 
higher resolution. The initial isotropic turbulent spectrum of fluctuations is chosen to be 
close to k""^ with random phases in both x and y directions. The choice of such (or even 
a flatter than -2) spectrum does not influence the dynamical evolution as the flnal state in 
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our simulations progresses towards fully developed turbulence. While the ISM turbulence 
code is evolved with time steps resolved self-consistently by the coupled fluid motions, the 
nonlinear interaction time scales associated with the plasma 1/k ■ Up(k) and the neutral 
1/k ■ V„(k) fluids can obviously be disparate. Accordingly, turbulent transport of energy in 
the plasma and the neutral ISM fluids takes place on distinctively separate time scales. 

Because of the different nonlinear time scales associated with the plasma and neutral 
fluids, mode structures can be different in the two fluids when they are evolved together and 
in isolation (i.e. decoupled). The former is shown in Figs ([3]) & Q- The flgures respectively 
show the evolution of various physical quantities in the plasma-neutral coupled system for 
a typic al set of ISM parameters a s descr i bed in the captions (for t y pical ISM parameters , 



seealso lZank I fll999h : IZank et al. I fll996[ ): iHeerikhuisen et al I (120061 ): iFlorinski et al I (12003 



20051 )). In the absence of charge exchange (see Fig. ([3]), right column), the plasma fluid 
evolves as an i deal MHD fluid and is known to typically develop sheet-like structures in the 
magnetic fleld (JBiskamp II2003I ). By contrast, the evolution of the magnetic fleld is modifled 
substantially when the two fluids are coupled through charge exchange. For instance, the 
sheet-like structures present in the turbulent magnetic field of the uncoupled plasma system 
instead are smeared in the coupled plasma-neutral system on the typical sheet-formation 
time-scales of ideal MHD. The neutral fluid, under the action of charge exchange terms, 
tends to modify the cascade rates by isotropizing the ISM turbulence on a slower time 
scale. This point is further elucidated in the subsequent section based on Kolmogorov-like 
phenomenological arguments. However if the coupled system is evolved further, sheets do 
form eventually in the magnetic fleld (after a relatively long time compared to the ideal 
MHD time-scales) and a turbulent equipartition is set up in the plasma and the neutral fluid 
modes. In any event, the small-scale sheet-like structures in magnetic fleld compress (or 
pinch) the plasma density. Hence the density fluctuations develop identical (thinner than) 
sheet-like structures that co-exist with small-scale turbulent fluctuations in their spectrum 
as shown in Fig. ([3]). The neutral fluid, on the other hand, evolves isotropically as stated 
above by forming relatively large-scale structures (see Fig. (jl])). 



5. Energy Spectra 

Spectral transfer in ISM turbulence progresses under the action of nonlinear interactions 
as well as charge exchange sources. Energy cascades amongst turbulent eddies of various 
scale sizes and between the plasma and the neutral fluids. In a freely decaying case, plasma 
and neutral fluids evolve under the influence of charge exchange forces which dramatically 
affect the energy cascades in the inertial range. This is evident from Figs ([5]) & ([6]) where 
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Kolmogorov-like fully developed turbulent spectra in the inertial range are shown respectively 
for the coupled neutral and plasma fluids. The neutral fluid energy spectrum in the inertial 
regime for the coupled system exhibits a fc~" spectrum, where the spectral index a ~ 3.2. 
On the other hand, the spectral index for plasma magnetic and kinetic energy spectra for the 
coupled system is a ~ 2.4. The inertial range spectral indices for the coupled plasma-neutral 
system in our simulations show a significant deviation from their corresponding uncoupled 
analogues which are respectively —3 and —1.7 for the neutral and plasma fluids, also plotted 
in Figs ([S]) & ([S])- The inertial range turbulent spectra are much steeper for the coupled ISM 
plasma-neutral system than the corresponding uncoupled or isolated neutral and plasma 
fluids in the decaying turbulence regime. Consequently, the steep turbulent spectra in the 
coupled ISM plasma-neutral system gives rise to an enhanced spectral transfer amongst 
the inertial range modes, and the energy cascade rates associated with the plasma and 
neutral (coupled) fluids can be understood by invoking a Kolmogorov-like phenomenology 
as described in the following. 

The typical nonlinear interaction time-scale in ordinary (i.e. uncoupled) plasma and/or 
neutral turbulence is given by 

Tnl ~ — ~ {kv^Y^, (7) 

where v^ or v^ is the velocity of turbulent eddies. In the presence of charge exchange in- 
teractions, the ordinary nonlinear interaction time-scale of fluid turbulence is modified by a 
factor fcce/fc such that the new nonlinear interaction time-scale of ISM turbulence is now 



k kvk 



The main justification of using the above expression originates from the fact that turbulent 
cascade is determined typically by the nonlinear interaction time (r„; ~ 1/kvk, where r„; is 
the eddy interaction time, k is characteristic turbulent wavenumber, and Vk is characteristic 
speed of eddy in wavenumber space) associated with the convection of turbulent eddies in 
the fluid momentum equation. Note that this time scale follows from Kolmogorov's phe- 
nomenology (Kolmogorov 1941) and corresponds essentially to the eddy interaction time in 
an ideal, non interacting, and non magnetized hydrodynamic-like turbulent fluid. Such time 
scales are essential to estimate spectral energy transfer (or cascades) rates {e ~ E{k)/Tni, 
where E{k) is turbulent energy per unit mode) in a hydrodynamic turbulent fluid. For a 
non-ideal, interacting and magnetized fluid, incorporating complex interactions amongst the 
turbulent eddies with a background or local mean magnetic field, the nonlinear eddy inter- 
action time needs to be modified. This was first pointed out by Kraichnan (1965) in the 
context of magnetohydrodynamic (MHD) fluids where the interaction of turbulent eddies 
with Alfven waves, excited as a result of a mean magnetic field, was addressed. Owing to 
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the presence of waves in a MHD fluid, turbulent correlations between velocity and magnetic 
field and the corresponding energy transfer time are determined primarily by r ~ (bok)^^, 
where 60 is a typical amplitude of a local magnetic field and is dimensionally identical to 
that of the velocity field by virtue of Elsasser's symmetry (Kraichnan 1965). Note carefully 
the modification in the energy transfer (or cascade) time scale because of the wave-turbulent 
eddy interaction process. Since the convective time is the only nonlinear time scale involved 
in the spectral transfer of turbulent energy, it it this time that accounts entirely for the 
nonlinear cascade in the coupled plasma-neutral heliosheath turbulence. Likewise, we in- 
troduce an analogous modification in the coupled plasma-neutral ISM turbulence where the 
energy cascade time scale due to charge exchange interactions amongst turbulent eddies, is 
determined explicitly by the factor kce/k (where the symbols have their usual meanings). 
It is this factor that determines how plasma and neutral fluids are coupled in heliosheath 
turbulence. Accordingly, the nonlinear spectral transfer time or the energy cascade rate is 
modified and is proportional to the factor kce/k such that the new nonlinear eddy interac- 
tion time in the coupled plasma-neutral heliosheath turbulence is t^l ~ {kce/k){l/kvk). The 
interaction time ttvl obtained in this manner for our coupled plasma-neutral problem not 
only adequately reproduces all three different regimes of interactions, but it also consistently 
describes the underlying physics in the context of heliosheath turbulence. 

On using the fact that kce is typically larger than k (or kc, the ch aracteristic turbulent 



mode , as defined elsewhere in the paper), i.e. k^e/k > 1 in the ISM (JFlorinski et al 1 12003 



20051 : IZank I Il999l ). the new nonlinear time is kce/k times bigger than the old nonlinear 
time i.e. t^l ~ {kce/k)Tni. This enhanced nonlinear interaction time in ISM turbulence is 
likely to prolong turbulent energy cascade rates. It is because of this enhanced or prolonged 
interaction time that a relatively large spectral transfer of turbulent modes tends to steepen 
the inertial range turbulent spectra in both plasma and neutral fluids. By extending the 
above phenomenological analysis, one can deduce exact (analytic) spectral indices of the 
inertial range decaying turbulent spectra, as follows. The new nonlinear interaction time- 
scale of ISM turbulence can be rearranged as 

kceVk 1 T^i ,„x 

kVk kVk Tee 

where Tee ~ {kceVk)~^ represents the charge exchange time scale. The energy dissipation rate 
associated with the coupled ISM plasma-neutral system can be determined from e ~ Ek/rNi, 
which leads to 

e^^^^^. (10) 

Keel ri Vk Kce 

According to the Kolmogorov theory, the spectral cascades are local in fc-space and the 
inertial range energy spectrum depends upon the energy dissipation rates and the character- 
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istic turbulent modes, such that Ek ~ e'^k'^. Upon substitution of the above quantities and 
equating the power of identical bases, one obtains 

Ek ~ ^2/3^-7/3 (^^) 

for the plasma spectrum (the forward cascade inertial range). The spectral index associated 
with this spectrum, i.e. 7/3 ~ 2.33, is consistent with the plasma spectrum observed in the 
(coupled ISM plasma- neutral) simulations (see Figs ([5]) & ([6])). Similar arguments in the 
context of neutral fluids, when coupled with the plasma fluid in ISM, lead to the energy 
dissipation rates 

This further yields the forward cascade (neutral) energy spectrum 

Ek ~ 5=^/3^-11/3 (^3) 

which is close to the simulation result shown in Fig. ([6]). 

A key issue finally is to understand what role charge exchange modes play in coupled 
plasma-neutral ISM turbulence as they result essentially from the nonlinear charge exchange 
interactions. To address this issue, we plot charge exchange sources associated with the 
momentum and energy equations, Eqs. ([I]) & ([2]), in Fig. ([7]). It appears from Fig. ([7]) 
that spectral energy is transferred predominantly at the larger scales by means of charge 
exchange mode coupling processes. The latter couples the large-scales, or smaller than kce 
modes, efficiently to low-fc ISM turbulent modes. It is primarily because of this kce — k mode- 
coupling in the smaller k part of the spectrum of the coupled plasma-neutral ISM turbulence, 
that energy is pumped efficiently at the lower k inertial range turbulent modes. The efficient 
coupling of the Fourier modes at low /c's further enhances the nonlinear eddy time-scales 
associated with the coupled plasma-neutral turbulence system which is consistent with the 
scaling Tjyi ~ {kce/k)Tnu where kce/k > 1, as described above. This consequently leads to 
the steepening of the inertial range spectra observed in Figs ([5]) & ([6]). By contrast, higher 
k modes, far from the energy cascade inertial range, are notably inefficient in transferring 
energy and momentum via charge exchange mode coupling interactions and are damped by 
small-scale dissipative processes in the coupled plasma-neutral ISM turbulence. 



6. Conclusions 

In conclusion, we have developed a self-consistent fluid model to describe nonlinear 
turbulent processes in a partially ionized and magnetized ISM gas. The charge exchange 
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interactions couple the ISM plasma and the neutral fluids by exciting a characteristic charge 
exchange coupling mode k^e, which is different from the characteristic turbulent mode k 
of the coupled system. One of the most important points to emerge from our studies is 
that charge exchange modes modify the ISM turbulence cascades dramatically by enhanc- 
ing nonlinear interaction time-scales on large scales. Thus on scales i > ice, the coupled 
plasma system evolves differently than the uncoupled system where large-scale turbulent 
fluctuations are strongly correlated with charge-exchange modes and they efficiently behave 
as driven (by charge exchange) energy containing modes of ISM turbulence. By contrast, 
small scale turbulent fluctuations are unaffected by charge exchange modes which evolve 
like the uncoupled system as the latter becomes less important near the larger k part of the 
ISM turbulent spectrum (see also Fig. ([7])). The neutral fluid, under the action of charge 
exchange, tends to enhance the cascade rates by isotropizing the ISM turbulence on a rel- 
atively long time scale. This tends to modify the characteristics of ISM turbulence which 
can be signiflcantly different from the Kolmogorov phenomenology of fully developed tur- 
bulence. We believe that, it is because of this enhanced nonlinear eddy interaction time, 
that a large spectral transfer of turbulent energy tends to smear the current sheets in the 
magnetic fleld fluctuations and further cascade energy to the lower Fourier modes in the 
inertial range turbulent spectra. Consequently it leads to a steeper power spectrum. It is to 
be noted that the present model does not consider an external driving mechanism, hence the 
turbulence is freely decaying. Driven turbulence, such as due to large scale external forcing, 
e.g. supernova explosion, may force turbulence at larger scales. This can modify the cascade 
dynamics in a manner usually described by dual cascade process. 

There are issues that need further exploration. For example, neutrals are observed to 
evolve isotropically in our coupled plasma-neutral ISM turbulence model. A relevant ques- 
tion, whether their coupling to the plasma fluid at large-scales reduces spectral anisotropy 
in the plasma fluctuations, remains to be investigated. Furthermore, our model is currently 
2D. It is then intriguing whether variations in the third dimension, i.e. 3D, introduce non 
trivial effects to our present 2D studies because turbulence in 2D differs considerably from 
that in 3D, particularly in the context of 3D neutral hydrodynamics fluid, where vortex 
stretching effect (which is absent in 2D) tears apart the large-scale structures and forbids 
the inverse cascade phenomena. By virtue of latter, 2D and 3D neutral fluids possess distinct 
spectral features characterized essentially by the number of the inviscid quadratic invariants 



(JBiskamp 1 120031 ) . How coupling of the plasma and neutral fluids in partially ionized ISM 



turbulence is modifled by the 3D interactions remains a subject of our future investigations. 

We add flnally that our self-consistent model can be useful in studying turbulent dynam- 
ics of partially ionized plasma in the magnetosphere of Saturn and Jupiter where outgassing 
from moons and lo and Bucephalus introduces a neutral gas into the plasma. 
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Heliosphere: region filled with solar wind plasma and bounded by HP. 
LISM/cloud: region beyond HP. 



HT: Heliotail 

HP: Heliopause 

IHS: Inner Heliosheath 

OHS: Outer Heliosheath 
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Fig. 1. — Schematic illustrating the plasma neutral simul a tion r egions both in the local 
interstellar medium and in the outer heliosheath (see IZank I (jl999l ) for a review of the solar 
wind interaction with the local ISM). 
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Fig. 2. — Schematic illustrating the significance of the charge exchange mode, which cor- 
responds to charge exchange interactions per unit area in the computational domain. This 
mode, being different from the characteristic turbulent mode k, is excited naturally when 
plasma and neutral fluids are coupled in the ISM by charge exchange. 
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Fig. 3. — (Left column) Mode structures of plasma when coupled with the neu- 
tral ISM fluid. The local ISM case is considered where typical LISM parameters are; 
7 = 5/3, nop,^^^ = 0.06cm-3, no„,„,,,i = 0.18cm-^,T = 6500K,Vth = I0km/s,a = 
6.9 X 10~^^, /o = 200AU, L(boxlength) = 5/o- (Right column) The corresponding quan- 
tities for an uncoupled ISM plasma simulation measured at a similar time step (i.e. t = 20 
in normalized units). Notice the magnetic field structures in the uncoupled case showing 
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Fig. 4. — The corresponding mode structures of neutral fluid in the presence of a plasma 
fluid. Simulation parameters are described in the previous figure. 
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Fig. 5. — The inertial range turbulent spectrum for a neutral fluid shows a fc~^'^ spectrum 
when coupled through charge exchange to a plasma fluid. The spectral index exhibited by 
this spectrum is larger than that of an ordinary (hydrodynamic) fluid. The latter exhibits 
a Kolmogorov-like k~^ spectrum in the forward cascade regime (drawn schematically). The 
numerical resolution for our spectral studies is 1024^ modes in 2D. 
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Fig. 6. — A similar steepening of turbulent spectra occurs in the plasma fluid when coupled 
to the neutrals through charge exchange. An uncoupled magnetofiuid plasma exhibits a 
Kolmogorov-like fc~^/^ spectrum in the forward cascade regime (drawn schematically). 
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Fig. 7. — Charge exchange spectrum associated with monientum and energy in both plasma 
and neutral fluids. Clearly the energy cascades are dominated by large-scale interactions, 
whereas small scales are not efficient enough to influence the turbulent spectra in the local 
ISM. 



